Search for: All records

Peptide coassembly offers novel opportunities for designing advanced nanomaterials. This study used coarse-grained molecular dynamics simulations to examine the coassembly of charge-complementary peptides, assessing various ratios and the role of charge and hydrophobicity in their aggregation. We discovered that peptide length, charge, and hydrophobicity significantly influence coassembly behavior, with more hydrophobic peptides exhibiting greater aggregation despite electrostatic repulsion. Beyond the coassembly of two peptides, we also observed that the coassembly of more than two peptides will likely lead to new assembly structures and properties. Our findings underscore the importance of peptide composition and length in tuning the coassembly and the resulting properties, thus facilitating the design of complex peptide nanoparticles for biomedical and biotechnological applications. 

The pKa values and associated protonation states of ionizable lipids in lipid nanoparticle (LNP) formulations are strongly dependent on their chemical environment. This phenomenon leads to poorly understood structure-function relationships that influence payload delivery, tissue-selective biodistribution, and manufacturing. For example, the charge- and biodistribution of an mRNA-loaded LNP can vary based on the type of ionizable lipid used, the molar ratio of its components, and its cargo. Yet, the spatial resolution of experimental protonation state measurements is currently limited to the apparent charge of an ionizable lipid averaged over all environments/conformations of an LNP — best represented by its apparent pKa value. Such measurements are too coarse to capture the heterogenous charge distributions of ionizable lipids in LNPs, which influence biocorona formation and interactions with the payload. Similar limitations are inherent to classical fixed protonation-state in silico models that cannot capture the environment-dependent protonation states and pKa values determining local pKa. To address this gap in experimental and computational tools available to accurately determine the local charge distributions in LNPs, this work now incorporates a scalable continuous constant pH molecular dynamics (CpHMD) model to simulate the self-assembly processes of five reported distinct LNP formulations. Parameters for ionizable lipids were generated from fitting fixed lambda-state calculations performed with Hamiltonian replica exchange (HREX) to improve conformational sampling during parameterization. Simulated systems were composed of 100 ionizable lipids (50 mol%), cholesterol (40 mol%), distearoylphosphatidylcholine (10 mol%), and mRNA (20 nucleotides) to model the interior of an LNP. Self-assembly was simulated for 100 ns at different pH values to validate the apparent pKa for each system. The theoretically calculated apparent pKa values matched reasonably well with those measured experimentally (mean absolute error = 0.5 pKa units), and all systems exhibited pH-dependent structures. Overall, this work provides a new computational platform technology to (i) predict the pKa values of ionizable lipids in different chemical environments and (ii) enable a structure-based way to model the heterogeneous, environment-dependent charge distributions of ionizable lipids in LNP systems typically encountered during LNP manufacturing and delivery. 

Wireframe DNA nanocages, an important type of DNA nanomaterials, exhibit exceptional programmability for chemical modifications, along with tunable size and shape. Nevertheless, the impact of their conformational fluctuations on cage design has not been thoroughly explored, despite speculation regarding its influence on biomedical applications. This study marks the first systematic examination of the conformational dynamics of prismatic DNA nanocages through molecular modeling and simulation. By comparing four different DNA nanocage topologies, we uncover design parameter combinations and conditions that facilitate access to varying conformational states. We observe the expansion and contraction of these cages across various topologies, hybridization states, and ionic environments (Mg2+/Na+), with their volumes varying from 15% to 150% of the ideal cage volumes. Our results indicate that the dynamics of DNA cages is influenced by the concentrations of Mg2+ and Na+ ions. Additionally, the flexibility of specific DNA strands can be manipulated, thereby altering the cage volume, through the selective hybridization of the cage edges. Ultimately, the conformational dynamics of DNA nanocages are captured in atomic detail. This study offers valuable modeling tools and methodologies to assist future DNA nanocage design endeavors. 

The passive membrane permeation of small-molecule drugs and relatively small hydrophobic peptides is relatively well understood. In contrast, how long polar peptides can directly pass through a membrane has remained a mystery. This process can be achieved with transcellular permeation enhancers, contributing significantly to the oral transcellular absorption of important peptide drugs like semaglutide — the active component in Ozempic, which is used as Rybelsus in a successful oral formulation. Here we now provide, for the first time, a detailed, plausible molecular mechanism of how such a polar peptide can realistically pass through a membrane paired with the permeation enhancer salcaprozate sodium (SNAC). We provide not only simulation results, obtained with scalable continuous constant pH molecular dynamics (CpHMD) simulations, but also experimental evidence (NMR, DOSY, and DLS) to support this unique passive permeation mechanism. Our computational and experimental evidence points toward the formation of permeation-enhancer-filled, fluid membrane defects, in which the polar peptide can be submerged in a process analogous to sinking in quicksand. 

Geometric isomerism in mechanically interlocked systems — which arises when the axle of a mechanically interlocked molecule is oriented, and the macrocyclic component is facially dissymmetric — can provide enhanced functionality for directional transport and polymerization catalysis. We now introduce a kinetically controlled strategy to control geometric isomerism in [2]rotaxanes. Our synthesis provides the major geometric isomer with high selectivity, broadening synthetic access to such interlocked structures. Starting from a readily accessible [2]rotaxane with a symmetrical axle, one of the two stoppers is activated selectively for stopper exchange by the substituents on the ring component. High selectivities are achieved in these reactions, based on coupling the selective formation reactions leading to the major products with inversely selective depletion reactions for the minor products. Specifically, in our reaction system, the desired (major) product forms faster in the first step, while the undesired (minor) product subsequently reacts away faster in the second step. Quantitative 1H NMR data, fit to a detailed kinetic model, demonstrates that this effect (which is conceptually closely related to minor enantiomer recycling and related processes) can significantly improve the intrinsic selectivity of the reactions. Our results serve as proof of principle for how multiple selective reaction steps can work together to enhance the stereoselectivity of synthetic processes forming complex mechanically interlocked molecules. 

Covalently linked molecular cages can provide significant advantages (including, but not limited to enhanced thermal and chemical stability) over metal-linked coordination cages. Yet, while large coordination cages can now be created routinely, it is still challenging to create chemically robust, covalently linked molecular cages with large internal cavities. This fundamental challenge has made it difficult, for example, to introduce endohedral functional groups into covalent cages to enhance their practical utility (e.g., for selective guest recognition or catalysis), since the cavities would have simply been filled up with such endohedral functional groups in most cases. Here we now report the synthesis of some of the largest known covalently linked molecular tetrahedra. Our new covalent cages all contain 12 peripheral functional groups, which keep them soluble. They are formed from a common vertex, which aligns the hydrazide functions required for the hydrazone linkages with atropisomerism. While we previously reported this vertex as a building block for the smallest member of our hydrazone-linked tetrahedra, our original synthesis was not feasible to be carried out on the larger scales required to successfully access the larger tetrahedra. To overcome this synthetic challenge, we now present an improved synthesis of our vertex, which only requires a single chromatographic step (compared to 3 chromatographic purification steps, which were needed for the initial synthesis). Our new synthetic route enabled us to create a whole family of molecular cages with increasing size (all linked with hydrolytically stable hydrazone bonds), with our largest covalent cage featuring p-quarterphenyl linkers and the ability to encapsulate a hypothetical sphere of approximately 3 nm in diameter. These results now open up the possibility to introduce functional groups required for selective recognition and catalysis into chemically robust covalent cages (without blocking the cavities of the covalent cages). 

The lack of biologically relevant protein structures can hinder rational design of small molecules to target G protein-coupled receptors (GPCRs). While ensemble docking using multiple models of the protein target is a promising technique for structure-based drug discovery, model clustering and selection still need further investigations to achieve both high accuracy and efficiency. In this work, we have developed an original ensemble docking approach, which identifies the most relevant conformations based on the essential dynamics of the protein pocket. This approach is applied to the study of small-molecule antagonists for the PAC1 receptor, a class B GPCR and a regulator of stress. As few as four representative PAC1 models are selected from simulations of a homology model and then used to screen three million compounds from the ZINC database and 23 experimentally validated compounds for PAC1 targeting. Our essential dynamics ensemble docking (EDED) approach can effectively reduce the number of false negatives in virtual screening and improve the accuracy to seek potent compounds. Given the cost and difficulties to determine membrane protein structures for all the relevant states, our methodology can be useful for future discovery of small molecules to target more other GPCRs, either with or without experimental structures.